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ABSTRACT 

The renormalization effects from the b-quark scale down to the non-perturbative 
QCD regime are studied for rare -B-decays at the heavy b-quark limit. Phe- 
nomenological consequences of these effects are investigated. We find that the 
anomalous scaling behavior plays a positive role in making non-perturbative 
model calculations consistent with recent CLEO measurements of B — > K*j. 



•k address after 1st Oct., 1993: Department of Physics, University of Tasmania, Hobart, 
Australia 7001 



1 Introduction 



Rare -B-decays are vital testing grounds for the Standard Model and therefore have 



received a lot of theoretical attention B y, ^7|. Besides examining the electroweak theory 



at the one-loop level, and providing quantitative information on the yet undetermined 
top quark mass and CKM matrix elements V t d,V ts ,V t b, radiative rare -B-decays could 
also be sensitive to new physics beyond the Standard Model. Recently, the CLEO group 
has reported the first observation of exclusive decay mode, B —>■ K* / y, and the measured 
average branching ratio, BR(B — > K*^) = (4.5±1.5±0.9) x 10 -5 , as well as the inclusive 
upper limit BR(B — > X s , -f) < 5.4 x 10~ 4 , are consistent with the Standard Model |L|]. 

To present theoretical estimates for exclusive processes, hadronic matrix elements 
governed by nonperturbative physics have to be evaluated. A number of effective ap- 
proaches on the basis of symmetry considerations and phenomenological models for QCD 
in the nonperturbative regime have been developed and applied to calculations of ma- 
trix elements for rare -B-decays @, ||, || |7|, [8|, || [10], [H| . As is well known, certain of 
these methods, for example, the nonrelativistic quark model, cannot take into account 
correctly the strong interaction effects at scales much larger than the QCD scale Aqcd- 
For current operators containing the heavy quark, the loop corrections with gluonic mo- 
menta between Aq C d and tuq lead to large logarithms of the type a s In A ™^ o ■ With 
renormalization group techniques, these large logarithms can be summed to all orders, 
resulting in an anomalous scaling factor. This kind of scaling factor has been worked 
out for the cases such as heavy meson decay constants Jm, B° — B° mixing|Tr2|, heavy 
meson semileptonic decay form factors [13|], and hadronic -B-decays [Q. In this letter, 
we will consider the anomalous scaling behavior from the b-quark scale down to the 
non-perturbative regime for current operators appearing in rare -B-decays upto next-to- 
leading order at the heavy b-quark limit and present the consequences of this scaling for 
exclusive rare -B-decays. 

In sect. II, we briefly review the renormalization from Mw to nib for rare decay 
operators in the full QCD with five quarks . More or less, we follow a similar procedure 
in discussing the anomalous scale below nib, which is presented in sect. III. The conse- 
quences of this QCD scaling for exclusive rare radiative -B-decays are shown in sect. IV. 
We summarize in sect. V. 



2 Renormalization in the full QCD 

To illustrate the general procedure of analyzing the anomalous scaling behavior through 
the renormalization group techniques, we briefly review the renormalization of current 
operators from heavy particle scale, such as the top quark and W boson, down to the 
bottom quark scale in QCD with five flavors. At the moment, we ignore the running of 
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the strong coupling constant between m t and Mw For the sake of concentration, we 
consider the rare radiative bottom quark decay. With no strong interaction in the form 
of QCD, the effective hamiltonian for the rare radiative bottom quark decay, after the 
t-quark and W-boson are integrated out from the (1-loop) photonic pengiun diagrams, 
reads 



n 



eff 



V2 



{s 3 + s 2 e* d )A(x t )(- 



with x t = m^/M^r. Here the contributions from virtual up and charm quarks have been 
ignored and a term proportional to the strange quark mass is dropped considering that 
mi, ^> m s . A(x) is the Inami-Lim function with the form 



r 2^2 



A(x) = x 



rp^ __|_ ^ ry> '_ ( O" T*l 1 Tl 1* 



[x 



(x-iy 
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The QCD effects play an important role in enhancing the radiative rare -B-decay, 
where the usual GIM suppresion factor in FCNC processes, is replaced by a logarith- 



mic one ln(^). Switching on strong interactions, the effective hamiltonian at <C M\ 



is determined by a proper operator basis. Following the notation of ref. ||, we have 



U 



AG 



eff 



v2 7= i 
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where 



o 1 = 


(s La ^b La )(c Lf 3'y^c L [3), 


o 2 = 
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(s La ^b La ) EqiqLpl^qLp), 
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o 6 = 




7 = 
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{^t)SL^T^b Rp Gl v . 
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with a and (3 the color indices. 

As the effective hamiltonian is renormalization scale independent, changes of the 
coefficients with the scale should be compensated by changes of these operators. This 
leads us to a renormalization group equation 



d 8 

a f jj i=i 



(5) 
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where 7 is the anomalous dimension matrix of the operator basis {Oj}. Introducing the 
QCD /9-function in terms of 



the general solution of the above equation can be written in a matrix form as 

C(M w /iM,g(jM)) = e X p[f BW dg^l\C(l,g(M w )). (7) 

Jg(M w ) fj{g) 

Explicit solutions can be obtained perturbatively. In the leading logarithmic approx- 
imation in which the terms like (a s InMjy/ n) n are summed to all orders, the /3-function 
and anomalous dimensions calculated at the one- loop level of QCD are needed, while 
matching conditions are of the zeroth order. For the four-quark operators 0\ — Oq, the 
matching conditions determined by the tree-level W-boson exchange are 

Cj(M w )=0 (j = 1,3,4,5,6), (8) 

C 2 (M W ) = 1. (9) 
The values of C 7 {Mw) and C%(Mw) follow from (one-loop) penguin diagrams. They are 

C 7 (M W ) = - l -A{x t ), (10) 

and 

C S (M W ) = ~D{x t ), (11) 



with 

D(x) = - 
v ; 2 



\x 2 - \x - 1 3a:lna: 



(x-if (x 111 ■ (l2) 

Furthermore at the next-to-leading order, the /9-function and anomalous dimensions are 
calculated to two loops of QCD and matching conditions are of the first order. In this 
case, the logarithmic terms like a s (a s \nM w / fi) 71 are summed to all orders. 

For the two-body decay B — > X s ^, which is modelled by b — > S7 at the quark level, 
only the photon magnetic moment type operator 7 contributes and the QCD corrected 
inclusive rare radiative decay rate is as follows 

T(B - X sl ) = ^j^(s 2 2 + s 2 3 + 2s 2 s 3 c 5 )\C 7 (ri\ 2 (13) 
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Practically, this rate may be normalized to the semileptonic 5-decay. In this way we 
remove the quark mixing angles in the small mixing limit and reduce the dependence 
on the bottom quark mass. So we have 



T(B -> X sl ) 



6a \CM\ 2 
% f{m c /m b ) 



1 



2a s (m b ) 
3tt 



-i 



g(m c /m b ) T(B -> X c lv{) 



(14) 



where f(m c /m b ) = 0.45 and g(m c /m b ) = 2.4 correspond to the phase space factor and 
the one-loop QCD corrections to the semileptonic decay, respectively. 

In the leading logarithmic approximation and with an anomalous dimension matrix 
for the truncated basis Oi, 2 , and 7 



7 



167T 2 



1 3 







3 -1 232/81 
16/3 



(15) 



the renormalization coefficient for the operator Oj has the following analytic form|| 
C 7 (fi) = -r/ 6 / 23 S^A(x t ) + [(58/135)( V - 10/23 ~ 1) + (29/189) (rT 28/23 - 1)]} , (16) 

with rj = a s (Mw) / ots(li) ■ Ref. has presented a discussion on the validity of the extra 
approximation with the truncated anomalous dimension matrix and estimated that the 
error in this coefficient function is less than 15%. Furthermore the calculation of ref. 



15[ in the leading logarithmic approximation without the truncation of the anomalous 
dimesion matrix also does not manifest significant change in the above estimates. Re- 
cently, the next-to-leading logarithmic effects which weaken the QCD-corrections to the 
b — > sj, have been included partly in ref. fl6| and result in a change of 10-20% in the 
inclusive decay rate. 

We notice that the upper bound of theoretical estimate for the inclusive rate is rel- 
evant to the renormalization scale which is usually chosen about the b-quark mass. 
Actually, based on eq. (14) and (16), a renormalization scale extrapolated almost to 
2.0GeV is used to set the upper bound of the estimate []1T 



BR(B 

in the minimal standard model. 



X s + 7) = (2 ~ 5) x 10 



— i 



(17) 
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3 Anomalous scaling at the heavy b-quark limit 



In this section, we present an analysis for renormalization of operators considered in the 
previous section below the bottom quark mass. There are three flavors of quarks in the 
area concerned. The masses of bottom and strange quarks lie somehow on the upper and 
lower bounds, respectively, whilst the charm quark locates in the middle. Based on this 
observation, we will scale down by treating the bottom quark as a heavy particle and 
integrating it in using an effective theory with four flavors. However, when we proceed 
into the regime between m c and m s , where both bottom and charm quarks are heavy, 
the light flavor number becomes three. 

Working in an effective theory where the bottom quark is treated as a heavy particle, 
we have the following expansion for Oj 

Oj = Y,D 3 i{m b /^ g(iM))Oi(n) + 0{—), (18) 
i mh 

where {0[} is an operator basis of the lowest possible dimension. Contributions of higher 
dimension operators vanish at the heavy bottom quark limit. Operators 0[ only depend 
upon the light scale ft and the large logarithm In ^ is transfered into the cofficients Dj h 
which obey the following renormalization group equations 



0, (19) 



where apart from different light flavor counting, the QCD /3-function retains the same 
form as the scale goes down, namely 



M = -g 



(20) 



with the one- loop (two-loop) coefficient listed in Table 1. Along with 7', the anomalous 
dimensions of {0[} in the heavy quark effective field theory (HqEFT), there is also a term 
in the above equation associated with 7, the amomalous dimensions of Oj in the higher 
scale theory. As shown explicitly in the appendix, the combination of the coefficients Cj 
in eq. (7) with the part of Dji associated to the anomalous dimensions 7^ gives one the 
anomalous scaling from M\y down to m&. We will drop the 7-term in eq. (19) when we 
work below the bottom quark mass. 

Since the photon magnectic moment type operator O7 dominates the rare radiative 
decay b — > s 7, we concentrate on it at the present stage. Let us consider the proper 
operators of dimension three [] in the effective theory where, to remove the heavy mass, 

1 At the moment we omit the phton field in O7 and the bottom quark mass that serves as the 
matching scale and is going to be set by the solution of fh(mb) — m&, with m(fi) the running mass. 
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the bottom quark field is redefined by 



h v (x) 



(21) 



Here is the four- velocity of the heavy bottom quark, which is a conserved quantity with 
respect to the low-energy QCD. Obviously, the first operator should be 0[ = SLa^ v h V R. 
Constructing from the Dirac matrix 7^ and the bottom quark velocity which is an exter- 
nal parameter also gives us the antisymmetric current operator 0' 2 = SL( , y^v u — v tl j u )h V R. 
Because the large logarithmic term arises only from the vertex of the heavy quark with 
the gluon, which becomes ig\ a v^ in the heavy quark limit, the renormalization is inde- 
pendent of the spin structure of the inserted current operators and does not mix them 
up. In this case, the anomalous dimension matrix reduces to a constant one and the 
renormalization group equation for the coefficients has a quite simple form 



4t + /%')J,-7V) 



dm = 0. 



<9/i dg' 

Expressing the solution of this equation in the same manner as eq. (7) gives 



(22) 



D 



71 



exp 



9 '(m c ) jW jV) 

g' {mb ) y P'(g') W) 9 P"(g") 



D n (m b ). 



(23) 



The second part of this solution arises when we cross the charm quark mass below which 
the light flavors reduce to three. The charm quark is only involved as a virtual heavy 
particle in loops and except different QCD couplings and light flavor numbers being 
used, 7" is the same as 7' which has the form 



a' 2 a' 2 
7 =7o T6^ +7l( 16^ ) + 



(24) 



with the so-called "hybrid "anomalous dimension of ref. |L2] and the two-loop coefficient 
shown in Table 1. 

As a first step we work in the leading logarithmic approximation and get the solutions 



D7M 



~a' s (m b )~ 


6 

25 


a'l(m c )~ 


2 
9 


_a' s (m c )_ 




[ . 





(25) 



and 



D 7l (n) = 0, for 1^1. (26) 

Here the matching conditions Dn(mb) = 1, and Dji{m}j) = for / 7^ 1 have been used. 
As matrix elements are concerned, the factor containing /1 should be cancelled by the 
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/z-dependence of hadronic states which are relevant to nonperturbative effects. This 
leaves us the anomalous heavy quark mass dependence 



lL LLA 



a' s (m c ) 



_6_ 

" 2,-) 



[a"(m c )\ 



(27) 



Then we consider next-to-leading logarithmic contributions. Using the matching 
condition for the current sTb to the g 2 order O, HI 



[r 



9 2 {m b ) 
24vr 2 



7 A ^ h\ 



(28) 



m b 



we find that at \i 

<V(75) — ► <v(7b), (29) 
just as at the zeroth order up to a vertex renormalization factor. Thus we also have only 
one non-zero coefficient, i.e. -D71 at the next-to- leading order. Once again we expand 
this coefficient function at m b in the QCD coupling 

J2 



Dn(m b ) = 1 + dx——; + 

lOTT 



(30) 



(31) 



and find the next-to-leading correction to the heavy mass dependence 

lai + /? ) — k — h K — ; . 

47T 47T 47T 

Here a /^-dependence factor that is absorbed into the matrix element is understood. 
It is worthwhile to point out that the combination d\ + k' is renormalization-scheme 
independent. The K-parameters determined by 

(32) 

2p 7o Po 

are listed in Table 2. Matching the vertex renormalization to the one loop in the full 
theory ]Tj| with that in the effective theory [|13|, [19| 



(Z Ml - 1) = (Z eff - 1) + 
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/2 



16n 



(33) 



at /x = m b in the MS* scheme gives us d\ = —Qcf, with cf = 4/3. 

Finally let us combine eq.(31) with eq. (27) and present the anomalous heavy quark 
mass dependence for the photon magnetic moment type operator O7 with both of the 
leading and the next-to-leading logarithmic contributions 



Q 1 



a s {m b ) 
a' s {m c ) 



25 



[a" s {m c )\ 5 [1 - 0.710c^(m b ) + 0.073a' s (m c ) - 0.060a;'(m c )] . (34) 
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The QCD fine structure constant accurate to the second order is 



a s {m 2 ) 



Air 



/?oln(m 2 /A^) 



1 - 



Alnln(m 2 /Af^) 
/3 2 ln(mVA|^) 



(35) 



As the heavy quark mass is concerned, we use the running mass in the MS scheme 
and take tjiq as the solution of m(mg) = mq. In terms of the invariant mass rh we have 
the following form pi 



m(n) = m [ln(///A 



MSI 



" ml) 

' 2/3 



1 - 



7m0 A lnln(/i 2 /A. 



+ 



2/? /3 2 ln(/i 2 /A§^) /3oln(/i 2 /A^ 



A/5' 7 



(36) 



where 7 m o (7mi) is the one- loop (two- loop) coefficient of the mass anomalous dimension 
(see Table 1). The values of K m in terms of eq. (32) are given in Table 2. With invariant 



masses in ref. [21] and h-jja ec L ua l to 0.25 GeV, we obtain scale masses in the MS 



scheme, = 4.39 GeV and m c = 1.32 GeV. Fine structure constants in Q T are 
a f s {m b ) = 0.204, a' s {m c ) = 0.332, a"(m c ) = 0.300, 



which give the following value 



T 



1.12 x 1.31 x (1 -0.139) S 1.26 



(37) 



(38) 



Several comments are in order: 

* The heavy quark mass dependence in the leading logarithmic approximation is 
universal to the current operators sTb with T any matrix in the Dirac space. However, 
the T-dependence in the full theory enters into the effective theory at higher orders 
through matching conditions. 

* Next-to-leading corrections are twofold. One comes from two-loop anomalous di- 
mensions and /^-function and does not depend on T. Another arises from matching 
conditions at the one-loop level and impacts on the correction differently for different 
T. For instance, the counterpart of d% in eq. (30) for vector and axial current operators 
is — 3cp |2(], [BJ. The large value of d\ for the present case leads to the result that the 
next-to-leading corrections are dominated by the a' s {rn^) term in eq. (34) that weakenes 
the leading effects. 

-k As we know, the photon magnetic moment type operator does not mix into four 
quark operators and the gluon magnetic moment type operator. The QCD corrections 
for the case in hand still remain unchanged when the mixings of these operators with 
O-j are taken into account. But these mixings add extra contributions, which have not 
been touched in this work and are believed not to change the essential feature of the 
anomalous mass factor. 
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4 Phenomenological consequence 



When a non-perturbation method at the QCD scale is used to evaluate hadronic ma- 
trix elements, the heavy-quark-mass dependence arising from renormalization should be 
taken into account. A number of such methods, such as the constituent quark model 
(CQM), MIT Bag model, and the heavy quark limit (HQL), have been employed in 
calculating rare -B-decay matrix elements |], |8], ^ [TI], |TT|. As a phenomenological appli- 
cation of our result in the previous section, we modify these calculation by incorporating 
the anomalous scale below the bottom quark mass. Our improved estimates for the ex- 
clusive rare radiative -B-decays are listed in Table 3. The Rq and R' are the ratio of 
exclusive to inclusive decay widths, i.e. 

F(B -> K* 7 ) , . 

before and after the modification, respectively. The top quark mass dependence and 
the coefficient CV are removed in this ratio. In the pure CQM and MIT Bag model 
calculations, the values of i?o are about 0.05. This yields an exclusive branching ratio of 
Br(B — > K*j) = 2.1 x 10~ 5 , for a top quark mass of 160 GeV. It is small comparied to 
the mean value of recent CLEO measurements. However, the anomalous scale factor of 
LLA increases this ratio almost upto 0.11, but the next-to-leading correction weakens 
this enhancement and results in a value of 0.082 for the improved ratio, which gives 
a branching ratio of 3.9 x 10 -5 . We find that incorporating the anomalous scale mass 
factor plays a positive role in making CQM and bag model calcualtions consistent with 
the CLEO date. Using the range of 3.5 ~ 12.2 for R in ref. [|1(]] obtained through varing 



the oscillator strength (3 between (3k = 0.34 GeV and (3b = 0.41 GeV, we get a range 
of (1.7 ~ 5.7) x 10~ 5 for the branching ratio of B — > K* / y decays. Our numerical results 
are slightly different from that of ref. |H| because of a different top quark mass being 
used. This range is lifted upto (2.7 ~ 9.4) x 10~ 5 by the anomalous mass factor. The 
improved estimates are completely in agreement with preliminay CLEO obseravtions. 

Even though the strange quark is not particular heavy compared to the QCD sale, 
HQL also presents reasonable estimate for the rare radiative -B-decay. We hope improve- 
ments will be included along this approach in the future. As the authors of ref. || point 
out, the large R at the third row of Table 3 is due to the use of a nonrelativistic recoil 
momentum. From the theoretical point of view, this is quenstionable considering that 
the recoil in the B to K* decay is very large. Experimentally, estimates with this large 
R for rare radiative -B-decays are not favored. 
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5 Summary 



In this letter, we have made an analysis of the anomalous heavy quark mass dependence 
of the rare -B-decays by considering the anomalous scale of current operators below 
the bottom quark mass. Combining this anomalous mass factor evaluated upto next- 
to-leading order with calculations using non-perturbative models, such as CQM and 
MIT Bag Model, gives us exclusive rare radiative -B-decay rates which are in excellent 
agreement with the recent CLEO measurements. On the other hand we observe that 
HQL also works as a preliminary approximation to B — > K processes. 

When we work below the bottom quark mass in this letter, mixings of four quark 
operators with the photon magnectic moment type operator have not been touched. We 
expect that they do not change the essential feature of the anomalous mass factor. The 
investigation of these mixings and their phenomenological consequences is under way 
and will be presented elsewhere. 
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Appendix A 

In this appendix, we present the formal solution of the renormalization group equa- 
tions, eq. (5) for \i < M and eq. (19) for /j, < m in the main text. We will show that 
the renormalization below m is effectively irrelevant to the anomalous dimension matrix 
of operators Oj. 

Once we write down the effective hamiltonian in terms of an operator basis 

j 

at /i <C M, the coefficient ought to satisfy the renormalization group equation, in a 
matrix form, 

(/^- 7 T )C = 0, (A2) 

where 7 is the anomalous dimension matrix. Introducing a transformation matrix V in 
such a way that 

(v-yv)« = tA, (^-3) 

where 7j are eigenvalues of the transpose of 7, we may construct an alternate opera- 
tor basis though the linear combination O = V t O. The effective hamiltonian can be 
rewritten in the same structure 

The advantage of this basis is that operators are multiplicatively renormalized and co- 
efficients have solutions like 

C,=exp[f%|]A'M- (^-5) 

Jg(M) p 

In the region /i < m < M where the proper operator basis is { 0[ }, there is the 
operator expansion 

I 

and the corresponding renormalization group equation has the form 

(^+7)D,-E-/H D fc = 0. ( A7 ) 

with 7' the anomalous dimension matrix of 0\. Using the transform (A. 3) and similarly 

(W~W T W) kl = %5 kh (A.8) 
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with 7j the eigenvalues of the transpose of 7', the coefficients have the following solutions 



D 



>ji = exp 



7,- 

sH P 



</M ,7,'' 



Djiim). 



(A9) 



In terms of these coefficients, the effective hamiltonian becomes 

Cj(M)exp 



jl jl { 



g( m ) 7- 

g(M) 13 



g'M ,7,' 

g'{m) 13' 



(A10) 

It is remarkable that the combination of the first piece of Dji, which is associated with 
the anomalous dimensions in the effective theory above m, with Cj in eq (A. 5) cancels 
the /x-dependence and gives the anomalous scaling behavior from M down to m, which 
is represented by the first coefficient of the effective hamiltonian in eq. (A. 10). The 
second coefficient, representing the anomalous scale from m down to /i, is determined 
by the anomalous dimensions in the effective theory at /i < m. 
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Table 1. Perturbation coefficients for the /3-function and anomalous dimensions 



/3-function 


anomalous dimensions 
of the quark mass 


anomalous dimensions 
of sTb in the HqEFT 


A) 




7m0 


Tml 


7^ 


7i 


11 - |JV> 


102 - f JV/ 


8 


8(M _ S N ) 
3V 2 3 JV // 


-4 


-(f + |7T 2 - f iV,) 



Table 2. Summary of values for the ^-parameter 



iV / = 4 


JV} = 3 




K l 


re" 


"'m 


if - -0.909) 


ff§(= 4 -° 6 ) 


g - gr^ -0.753) 


f°-(^3.58) 



Table 3. Improved branching ratios for 5 — > K*^ 



Model 


ref. 


Ro(%) 


R'{%) 


Br(B -> K*7) x 10 5 


CQM(r) 


i 


4.5 


7.4 


3.5 


CQM & MIT Bag 


i 


6.0 


9.8 


4.6 


CMQ(nr) 


i 


21. 


34. 


16.1 


HQL+CQM 


m 


3.5 ~ 12.2 


5.7-20.0 


2.7 ~ 9.4 


HQL 


i ii i 


12. ~ 27. 


20. ~ 44. 


9.4 ~ 20.5 
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